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Abstract 

0^ 

CN ' We present an integrable Hamiltonian which describes the sinh-Gordon model 

on the half line coupled to a non-linear oscillator at the boundary. We explain 
^ , how we apply Sklyanin's formalism to a dynamical reflection matrix to obtain 

\ this model. This method can be applied to couple other integrable field theories 

■ to dynamical systems at the boundary. We also show how to find the dynami- 
\Q . cal solution of the quantum reflection equation corresponding to our particular 

\ example. 

1 Introduction 

O . Integrable field theories in two dimensions provide us with a theoretical laboratory to 

study non-perturbative phenomena in high-energy physics. They also possess numer- 
ous applications in fluid dynamics, statistical physics, condensed matter physics, and 
^ ■ quantum optics. 

■ In recent years it has become possible to extend many of the techniques and results 
of integrable models to field theories with boundaries, i.e., theories defined either on 
the half line or on an interval, see for instance P|, EBl El, O, B] . 



We introduce a new class of integrable field theories with boundaries where, instead 
of imposing fixed boundary conditions, we couple the boundary field to a mechanical 
system. 

We begin by giving a very concrete example in section ^ where we describe the 
sinh-Gordon or sine-Gordon model coupled to a non-linear oscillator at the boundary. 
Then in section |] we explain more generally how one can couple an integrable field 
theory to a mechanical system in such a way that the integrability is not broken. We 
use a generalization of Sklyanin's technique for the construction of integrable boundary 
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conditions. The new ingredient is that the solution of the reflection equation is chosen 
to depend on boundary degrees of freedom. In section |^ we speciahze this technique 
to the case of the sine-Gordon theory, providing some details of how we arrived at 
the model described in section 0. Finally, in section ^ we describe how to obtain a 
particular dynamical solution of the quantum reflection equation. This is the solution 
whose classical limit went into the construction in the previous section. We end with 
discussions in section |[ 

Before we begin let us survey other works which deal with coupling to degrees of 
freedom at the boundaries. These fall into three categories: 

• Field theory: In H] the sine-Gordon model is coupled at the quantum level 



to a q-oscillator. In |^ it was observed that in order to derive the flxed boundary 
conditions for the supersymmetric sine-Gordon model from an action, it was convenient 
to introduce a fermionic variable at the boundary which could however be integrated 



out again immediately. In ||T^ a free fermion fleld theory is coupled to a dynamical 
boundary degree of freedom. 

• Integrable quantum spin chains: chains coupled to extra spins on the boundary 
are used in the study of Kondo impurities coupled to strongly correlated electron 
systems, see e.g. p8|. In 0, |lT| new dynamical solutions of the quantum reflection 



equation which are not of the RKR type were used to describe the coupling of the 
boundary spins. 

• Systems with a flnite number of degrees of freedom: Kuznetsov has used dynami- 
cal solutions of the classical reflection equation to couple mechanical tops to integrable 
non-linear lattices such as the Heisenberg chain or Toda lattices. Our work can be seen 
as the extension of these ideas to fleld theory while preserving integrability. 



2 The sinh-Gordon model coupled to an oscillator 
at the boundary 

In this section we present an integrable Hamiltonian describing the coupling of the sinh- 
Gordon fleld theory to a non-linear oscillator at the boundary. We leave the details of 
how we constructed this Hamiltonian for later sections. 

The sinh-Gordon model describes a relativistic 1+1 dimensional self-interacting 
massive bosonic fleld (f){x,t). The Hamiltonian of the sinh-Gordon model restricted to 
the half-line is 

/■° /I 1 - \ 

i^shG = J dx (^-vr^ + -(5.0)2 + -^(cosh/30 - 1)J . (1) 

Here vr is the conjugate momentum to the fleld 0, i.e., 

{7c{x),<P{x)} = 6{x-y). (2) 
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/3 is the real sinh-Gordon coupling constant and m sets the mass scale. We let space 
be the half line from x = — oo to x = 0. As usual the field (p is assumed to vanish at 
X = —oo but not at a; = 0. As will become clear from the construction in section ^ we 
could also have taken space to be an interval [x_ , x+] and placed a mechanical system 
at both ends. 

We describe the coupling of the sinh-Gordon field 0(0) at the boundary a: = to a 
non-linear oscillator through the boundary Hamiltonian 

Here q and p are the position and momentum variables of the oscillator and they 
obey the canonical Poisson bracket relation {p,q} = 1. The new free parameter M 
determines the mass of the oscillator. 

It is instructive to expand the Hamiltonian for the oscillator for small p and q. One 
obtains 

= 2W/ + ^^'^ + ^ ~* (^) + + • 

where 

= Me^'^(°)/2 and u = m/2 . (5) 

We see that the frequency u of the oscillator is fixed by the requirement of integrability 
to be equal to half the mass of the sinh-Gordon field. The effective mass of the oscillator 
depends on the value of the sinh-Gordon field at the boundary. The exact form of the 
higher non-linear terms in the Hamiltonian of the oscillator are fixed by integrability. 
To shorten the formulas we introduce the rescaled variables 

which have Poisson bracket {p, q} = /3^/4. In terms of these the boundary equations 
of motion are 

|p={ifosc,J>} = -ye^^(°)/^inhg, (7) 

|g = {//osc,g} = f e-^*(°)/^inhp. (8) 

To determine the equations of motion of the sinh-Gordon field (p one needs to use 
the full Hamiltonian 

H = HshG + Hose- (9) 



3 



One finds 



j^<P{x)={HA{x)} = Ti{x), (10) 
-n{x) ={H,nix)} = d'^x) - sinh/30(a;) (11) 

Note the term proportional to 6{x) in the equation of motion for tt{x). It has two 
sources: the dx4> arises when one performs a partial integration in the sinh-Gordon 
Hamiltonian which produces a boundary contribution at x = and the other two 
terms arise because the Hamiltonian of the oscillator contains 0(0). 

We are only going to allow solutions to these equations which are continuous on 
the left half-line [— oo, 0]. So in particular we require continuity of '7r(x) at x = 0. This 
implies that the S{x) term in the equation of motion for it{x) vanishes because this 
term would otherwise force Tr{x) to develop a discontinuity at x = 0. We therefore 
know that the solutions will satisfy 

9,0(0) = ? (cosh p e-^'^^o)/^ _ ^Qgi^ ^ g/3<A{0)/2 j _ (^2) 

We can combine the equations of motion ( p!OD and ([l^) to obtain the usual sinh- 
Gordon equation of motion 

2 

777 ^ 

9^0 -9^0 = -^ Sinh /?0. (13) 

This equation for together with the boundary condition (|T2p and the equations of 
motion (|^) and (P) for the oscillator form one system of coupled equations involving 
both ODE's and a PDE which needs to be solved to obtain the time evolution of the 
system. 

We call this system "integrable" because, as we will show in section |, we can 
construct an infinite number of conserved higher spin charges J„ which are in involution 
with each other. The first non-trivial charge beyond the Hamiltonian is the spin 3 
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charge 

(tt^ + 5{dAY) cosh/30+ ^ (cosh 2/30 - l) j rfx 
1 

— e'^'^^'^^^'^ - coshp cosh^ q coshp -tt^ cosh q 

\ 2 2 2m"' 

+ e-^'^*^")/^ Q coshp + i cosh^p^ 

/ 3 1 

_ g-/3<^(o)/2 -coshg cosh^pcoshg -vr^coshp 

\ 2 2 2m"' 

H — ^ sinh q sinh p . 

m 

The boundary condition (0) is similar in form to the previously known integrable 
boundary conditions O, IT! 



with the crucial difference of course that the parameters eo and ei were fixed numbers 
rather than dynamical variables as in our case. Only in one case do our boundary 
reduce to the fixed boundary conditions (p!5D, namely when the boundary oscillator is 
at rest at g = p = 0. This corresponds to the boundary conditions with eo = Ci = 1. 
The quantum fluctuations of the boundary oscillator will probably imply that our 
model never reduces to the usual boundary conditions in the quantum case. 

After quantization the sinh-Gordon model describes scalar massive particles. The 



direct solution of the quantum theory is rather difficult |25[. However because of the 
existence of higher spin conserved charges one knows that there is no particle production 
and that the particle scattering factorizes into a product of two-particle scattering 
processes. The corresponding scattering amplitude has been obtained by analytical 
continuation in the coupling constant from the breather scattering amplitude in the 
sine-Gordon model ||29l . 

To describe the sinh-Gordon particles on the half-line one also needs to give the 
reflection amplitudes. In the case of the fixed boundary conditions ([T5|) this amplitude 



can again be obtained from the corresponding breather reflection amplitude in the sine- 



Gordon model [|lj]. Because we expect that the results for our model can be obtained 



similarly we now turn our attention to the sine-Gordon model. 
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When we let the sinh-Gordon coupling constant (3 become purely imaginary, i.e., if 
we set (3 = i(3 with j3 purely real, then the sinh-Gordon Hamiltonian (|lD turns into the 
sine-Gordon Hamiltonian 

^sG = f_ dx Qtt^ + ^{d,<pf - ^(cos/30 - 1)^ . (16) 

Under the same replacement the Hamiltonian for the oscillator becomes 

Notice that in general this Hamiltonian is not real. We are familiar with this situation 
from imaginary coupled affine Toda theories. There the non-real Hamiltonian implies 
that the classical soliton solutions are complex. The energy of these configurations is 
nevertheless real [jl5|. Similarly here the classical boundary solutions will be complex 



but we expect that their energies will be real. 



3 Sklyanin's formalism 



In a seminal two-page paper p3[ Sklyanin described how one can impose boundary 
conditions on an integrable field theory without breaking integrability. Below we will 
describe how Sklyanin's formalism allows us to couple an integrable field theory to a 
mechanical system at the boundary rather than imposing a fixed boundary condition. 
This section will be easier to understand for readers who are familiar with the approach 



to integrable models described in |T0 



In the following we assume that there exists a pair of matrix valued functions 
ax{X, x) and at(A, x) which depend on the fields of the theory, their conjugate momenta, 
and on a spectral parameter A = G C, so that the classical equations of motion of 
the field theory are equivalent to the Lax pair equation 

[dx — a'xiX,x),dt — at{X,x)] = for all A. (18) 

Here ax{X, x) and at{X, x) depend on x and t only implicitly through their dependence 
on the fields, dx and dt denote total differentiation with respect to the space or time 
variable. If ax{X,x) and at{X,x) are thought of as the components of a connection 
then eq. (^) is the zero curvature condition for this connection. 

Equation ([T8| ) is the compatibility condition for the overdetermined system of equa- 



tions 



dT 

= ax{X,x+)T, (19) 



dx+ 
dT 

— = at{X,x+)T -Tat{X,x^ 
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where the transition matrix T = T{x+, x_, A) is defined to be a solution of the differ- 
ential equations (|19D with the initial conditions T(x_,a;„, A) = /. It can be expressed 
as the path ordered exponential of ax{X,x) from a;_ to Xj^ 

T(x+, A) = Pexp ( / aa;(A, (20) 

so that the operators at points nearer to x^ are further to the left. We assume that 
the Poisson brackets for the functions a^iX, x) can be written in the form 

12 12 

{a^{Xi,x),a^{\2,y)} = 5{x - y)[r{]n{Xi/\2)),a^{\i,x) + a^{\2,y)]. (21) 

1 2 

Here we used the short hand notations a(A, x) = a(A, x) ® /, a(A, x) = I a(A, x) and 
the r-matrix is independent of the field or its conjugate momentum. It follows that 

{T(x+, x_, Ai), T(x+, A2)} = [r(ln(Ai/A2)), T(x+, Xi)T{x+, A2)] . (22) 

For simplicity we assume below that the r-matrix has the property | 

r{e) = -r{-e). (23) 

Let us now introduce two matrix valued functions K±{6) of the spectral parameter 
A = e^, not depending on the fields. However, in departure from the situation described 
by Sklyanin in we let these boundary K matrices be dynamical. This means that 
we enlarge the phase space of the theory by introducing extra degrees of freedom which 
we think of as describing mechanical systems placed at the boundaries. We then let the 
K matrices depend on these new dynamical variables in such a way that their Poisson 
brackets are given by the following classical refiection equation: 

{K±iO),K±i9')} = [ri9-9'),K±ie)K±{9')] (24) 

+ k±ie)r{e + e')K±ie') - K±{e')r{e + e')k±{e). 

Furthermore we assume that the dynamical systems on the left and right boundary are 

1 2 

independent so that {K±{9), K^{9')} = 0. 

Following Sklyanin we define the functional 

r{x+,x^,X) = T{x+,x.,\)K^i\nX)T-\x+,x.,l/\) (25) 

which generalizes the transition matrix to the boundary case. Using eq. ( p3D and that 
the Poisson bracket is antisymmetric, satisfies the Jacobi identity and has the property 

^This could however be relaxed, see for example |^ for an application of Sklyanin's formalism to 
afRne Toda theory where the r-matrix does not possess this property. 
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{A, BC} = {A, B}C + B{A, C} it is straightforward to show that T(x+, x_, A) obeys 
a Poisson bracket relation similar to those for the K^s given in eq. p^). A family of 
transfer matrices is defined by 

r(A) =tr(K+(ln(A))T(x+,x_,A)). (26) 

These are in involution for any values of the spectral parameters 

{r(Ai),r(A2)} = for all (Ai,A2)gC (27) 

provided appropriate boundary conditions are imposed at x±. As the generating func- 
tion r(A) can be expanded about the singularities in the transition matrix T(x+, a;_, A) 
it gives an infinite number of quantities /„ in involution with each other. Among these 
we will identify the Hamiltonian of the model. Then it follows that the {/„} are time 
conserved. The model is thus integrable on the interval Different dynamical 

systems can be coupled at the two boundaries by choosing different K_ and Kj^. 



4 Derivation of the sine-Gordon example 

In this section we will provide some details needed to apply the general method de- 
scribed in the previous section to the sine-Gordon model in order to derive the coupling 
to a boundary oscillator described in section ^. 

The equation of motion (13) for the sine-Gordon field is representable as the zero- 
curvature condition for the Lax connection a^(A,a;) written in terms of the standard 
Pauli matrices ak, A; = 1, 2, 3, as 

a^{\x) = __a3 + -(A + -)sm(— )ai + -(A--)cos( — )a2, (28) 

= 4.5^"^ + 4^(^-a)^^"M^^ + 4^(^+a)^°'M^^- (2^) 

We define 7r(x, t) = ^^t^Jl, Notice that the matrix ax{X,x) possesses two singularities 
located at |A| = and |A| = oo. 

To extend the standard definition of the Poisson bracket, we change the range of 
the integrals from [— oo, +oo] to Then, using the notation 7r(x) = 7i{x,t) and 

(j){x) = 0(x, t) at fixed time t we define the Poisson bracket as 



x+ 



60i 6O2 6O1 60: 



dOi d02 dOi dO 



{0,,02}= I dx -^-^ --^-^ +^^~^^ (30) 
Jx ld7T[x) d(p[x) d(p[x) d7r[x)^ op oq oq op 

for any observable Oj. Obviously, this bracket possesses the basic properties of a 
Poisson bracket. It is skew symmetric and satisfies the Jacobi identity. Also, we 
have {Ci,C»2C3} = {Oi, 02)03 + 02{0i,0^}. Using the definition (|^), the non- 
vanishing Poisson brackets in the sine-Gordon field theory at constant time slices are 
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{n{x),(j){y)} = 6{x — y) and for the boundary variables p(t), q{t) we have {p, g} = 1. 
Now, by calculating the Poisson brackets of the entries of the matrix ax{\iX) it gives 
a unique solution to eq. (|2^) : 

/^^coshm. ^ N / ^ / ^ 

losmh(t^) losmh(t^) 



To construct the generating function ( P6[) of the integrals of motion for the sine-Gordon 
model with dynamical boundaries, we need (p, g)-dependent solutions K±{6) to the 
classical reflection equation (|2^) . Let us introduce the rescaled variables p = pP/2\^ 



and q = qP/2\^. Then the matrix 

^ (0) — 2 f ^^^Hp + ~ cosh(p - q)e~^ 2 sinh^(^) - 2 sinh^(p) \ , , 

\^ 2 sinh^(g) — 2 sinh^(^^) cosh(p — g)e^ — cosh(p + g)e^^y 

satisfies the classical reflection equation (^) with the classical r-matrix r{6) = rsai^) 
defined above. We could place additional degrees of freedom at the left boundary as 
well and use a of a similar form. However we would not learn anything new and 
for simplicity we choose {6) = I which is a trivial solution of the classical reflection 
equation (0) and then move the left boundary off to —oo, i.e., we restrict the sine- 
Gordon field theory to the half-line. We then assume that the field and its conjugate 
momentum satisfy the Schwartz boundary condition = and 7r(x_,t) = at 

X- = —oo. Using eq. (pSD the generating function (p6D becomes 



r(A) = tr(K+(ln(A))T(0,-cx),A)r(-oo,0,l/A)). (33) 

Fixing Xm{9) = |, the transition matrix T(— oo, 0, A) has singularities located at 
A = ±ioo and A = ±iO. There will exist two infinities of involutive integrals coming 
from the coefficients of the Laurent expansions about these two values. First, let us 
consider the asymptotic expansion of the transfer matrix r(A) as |A| oo. Substituting 
the solution K+{9) of (|32D into (|33D and expanding about A = ioo, we obtain a Laurent 
series in A which provides an infinite number of quantities /„: 

oo J 

ln(r(A)/A^)= El^- (34) 

n=— 1 

In performing this expansion it has helped us to look at how the corresponding calcula- 
tion was performed in [|l^] for the case of non-dynamical K. We find /_i = — ^^^|l^oo 
and Jo = 0. The next quantity Ji = —^H + const, gives the Hamiltonian H for the 
system, where 

H = l\x[e{-x){^[n' + id,<f^f]-'^icosm-l)') (35) 
+ 6{x)^ (cosh(p + q) e'^-^/' + cosh(p - q) e'^^^/^N 
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This can be seen to agree with the Hamiltonian given in section ^ after we perform the 
canonical transformation 



(36) 



The next term in the expansion (3^) of order 0(1/A^) reduces to /2 = after using 
the boundary condition (p!2D . Similarly the 0(1/A^) term in (plf ) reproduces the spin 
3 charge J3 given in (|I4|), again after the canonical transformation (pUf). It is believed 
that the existence of such higher local integrals of motion is a sufficient condition for 
the classical integrability of the system. 

Finally we observe that because the transfer matrix r(A) in (|33D is invariant under 
the simultaneous replacements A —l/\,(f){x,t) — > — 0(x, t) and q —q, the expan- 
sion of r(A) around the singularity at gives analogous quantities J„. Notice that H 
and I3 are invariant under this transformation. 

It is rather interesting to notice the following: If p and q are fixed c-numbers then, 
the (p, g)-dependent part of the Poisson bracket (^) disappears and the l.h.s. of eq. 
( p^ vanishes. Considering the two leading terms in the expansion in A of the classical 
refiection matrix K^{6) given in eq. (0) we obtain 



^+(ln(A)) 



I A|— >oo 



^2 cosh(p + g) 1 

—1 I cosh(p — q) 



0(A° 



(37) 



which can be shown to satisfy eq. (p^ . Let us now introduce the notation 



cosh(p + q) = P + iQ , 
cosh(p — q) = P — iQ . 



(38) 



Multiplying K_|_(ln(A)) by 1/A^ the reflection matrix now writes 



9 P HO 
ir+(ln(A)) ~ —1 + 2^2 + ^^3 + 0(1/ X^). 
A A 



(39) 



Up to order 1/A this is the same K matrix used in |T9[ to construct the integrable 



boundary conditions used by Ghoshal and Zamolodchikov in [|l^]. Of course, our next 
quantities I2 and /a differ from the ones in [|l^ due to the contributions of order 0(1/A^) 
that appear in our case. 
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5 Dynamical solutions of the quantum reflection 
equation 

In this section we look for operator valued solutions K,±{9\ a) of the quantum reflection 
equations 



R{e - e')t-{e + a/2; a)R{e + 6')k^{6' + a/2; a) 



h-{e' + a/2; a)R{e + 9')K-{e + a/2; a)R{e - 9') , (40) 



1 2 

R{-e + e')lC\{e - a/2; a)R{-e - e')JCX{e' - a/2; a) 

2 1 

= IC'4e' - a/2; a)i?(-^ - e')ICX{0 - a/2; a)i?(-^ + 6') . (41) 

Here R{0) is a given quantum i?-matrix. We see that the quantum reflection matrix 
)C-{6;a) for the left boundary has to satisfy a different equation to the quantum 
reflection matrix )C^{6;a) for the right boundary. The classical reflection equation 
can be obtained as a limiting case of either of the quantum reflection equations 
and (||T|)0. 



One useful way to think of the reflection equation (|40| ) is to view the entries of 
the matrix K,-{9; a) as the generators of an associative algebra with quadratic algebra 
relations given by (^0[). Such an algebra is often called a reflection equation algebra. 
Finding a solution of a reflection equation is therefore the same as finding a repre- 
sentation of the corresponding reflection equation algebra. Below in (|46| ) we will give 
an infinite dimensional representation in terms of position and momentum operators q 
and p. 

Operator valued reflection matrices have several applications: 1) they can be used 
in integrable lattice models or quantum spin chains to couple additional boundary spins 



to the model [11, 30], 2) they can describe the reflection of particle excitations from 
a boundary if there are degeneracies in the boundary spectrum, and 3) their classical 
limit can be used to construct integrable models p6| , p!7| . 

As an example, let us consider the following trigonometric and hyperbolic minimal 
solutions of the quantum Yang-Baxter equation: 

^(^,7) = ^(/®/ + a3®(T3) +^(/®/-a3®a3) 

+ ^ (CTI ®CTi+a2®(J2). (42) 



"Our /C± are related to the K± of |2| by JC±{e + a/2) = K±{e) and a = rj. 

^In the Hmit a ^ one recovers the classical r-matrix (up to an overall ^-dependent coefficient) 
R{9) ^ I + ar{9) + 0{a^) where / is the identity matrix. The classical reflection matrices if ± (6') are 
defined to be the first term in the expansion of the quantum refiection matrices: K,± (9] a) = K± (9) + 

a5K±{9) + 0{a^). 
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This i?-matrix gives the Boltzmann weights of the six vertex model which is known to 
be related to the XXZ and the XXX (in its rational limit) spin chains. It possesses 
three different regimes called antiferroelectric (I), trigonometric (II) and ferroelectric 
(III) with 

(I) a{e) = smh{^ - 6) b{e) = sinh{e), c(^) = sinh(7), 7 > > 0, 

(II) a{e) =sm{-f -6), b{e) = sm{e), c(^) = sin(7), vr > 7 > > 0, 
(III) a{e) = smh{e + 7), b{e) = sinh(^), c{e) = sinh(7), ^ > 0, 7 > 0, 

where the parameter 7 characterizes the anisotropy of the model. In particular, the 
trigonometric regime of the six- vertex model (II) describes the critical (zero gap) limit 
of the eight- vertex model 0. 

A dynamic solution of the quantum reflection equation corresponding to this R- 



matrix was already given in |16, 18 without derivation. We will derive a similar 



solution below but with some minor quantum adjustments. 

For further convenience, defining the K-mathx at the quantum level by 



\D{e) E{e)) 



(43) 



and using the expression for the i?-matrix written above, the reflection equation (0) 
now reduces to eight functional equations: 

(i) a_c+ {BD' - B'D) + a_a+ [A, A'] = 0, 



iiii 



h_h+{AE' - E'A) + c^c+[E,E'] + c^a+{DB' - D'B) = 0, 
a_c+ {DB' - D'B) + a_a+ [E, E'] = 0, 
(iv) c_b+ {EA' - E'A) + 6_c+ {AA' - E'E) + 6_a+ {BD' - D'B) = 0, 
(v) h-h+AD' + C-C+ED' + c^a+DA' - a^a+D' A - a^c+E'D = 0, 
(vi) h_a+BE' + c_h+EB' + b_c+AB' - a_b+E'B = 0, 

where we use the shorthand notations a_ = a{6 — 6'), = a{6 + 6' — a) and similarly 
for b and c as well as A = A{6) and A' = A{6') and similarly for B, D and E. The 
remaining two equations are obtained from (v) and (vi) through the substitutions 
A ^ E and B ^ D. Let us now focus on solutions of kind (III). To find the solution 
of the equations (i) — (vi), we assume the following form for the i^-matrix: 

B{e) = -2sinh2^ + 2f/ D{e) = 2smh^ 6 - 2V 

Setting 7 = a in case (III), the generators {F,G,U,V} have to satisfy the following 



12 



relations: 

[F, G] = -2 sinh a{U - V) , [f/, V] = ^^(f^ _ G') , 
FVe-'^ - VFe'^ - F sinh a + G sinh 2a/2 = , 

GVe"" - VGe''^ + G sinh a-F sinh 2a/2 = , (44) 
FUe-" - VFe" - F sinh a + G sinh 2a/2 = , 
Gf/e" - f/Ge~" + G sinh a-F sinh 2a/2 = . 

One can find a representation for the generators {F, G, U, V} satisfying these quadratic- 
linear relations in terms of the position and momentum operators q and p with com- 
mutation relation 

\p{t),q{t)]=a. (45) 

This leads to the following solution of the quantum refiection equation (^): 

, ._/cosh(p-g)e^+"/2_^,Qsh(p + g)e-^'-°/2 2 sinh2(g) - 2 sinh2(^) 

2sinh2(^) -2sinh2(p) cosh(p + g)e^+"/2 _ cosh(p - g)e-^^"/- 

(46) 

This /C_(6';a) satisfies (|40D also in regime (I) with 7 = —a. With the substitution 
6 —>■ i6 and 7 27 the solution corresponding to the regime (II) follows immediately. 
The matrix /C+(0; a) = lC'_{—9\ a) provides a solution of (^I]) . The classical refiection 
matrix (P2|) used in the previous section is obtained as the classical limit of /C_(6'; a). 



As was observed by Sklyanin, given a c-number valued solution K of the refiection 
equation one can always construct an infinite number of additional operator valued 



solutions by dressing K with /2-matrices to RKR, see p5[ for details. The solution 



which we derived above is probably not of this factorizable form. 



6 Discussion 

We have explicitly constructed an integrable classical field theory with extra degrees 
of freedom living at the boundary. We have shown how to obtain concrete expressions 
for the Hamiltonian and higher spin conserved charges. 

While in this paper the sinh-Gordon model coupled to an oscillator at the boundary 
was given only as an example of the kind of model one can obtain, we believe that it 
is of great interest in its own right and we intend to study both its classical solutions 
and its quantization. In both endeavors we will be helped by the large amount of work 
which has been done already on the sine-Gordon model with fixed boundary conditions, 
see for example f^^, |TB|, ^, ^, ITf. 

Of particular interest are the classical solutions describing oscillating boundary 
states. In the case of the fixed boundary conditions these so called boundary breather 
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solutions were found and quantized semiclassically to obtain the spectrum of boundary 
states 0, § . In our model we will look for classical solutions in which not only the field 
near the boundary but also the boundary variables p and q oscillate. It is possible 
that there will be two or more solutions with the same energy, they might for example 
be obtained from one another through the transformation (p —(t>iP ~^ P 

which is a symmetry of the theory. This question of how many degenerate states there 
are is very relevant to the determination of the reflection matrices. As in the case of 



fixed boundary conditions described in [^], the soliton reflection matrices will have to 



satisfy the quantum reflection equation. However if there is a degeneracy of boundary 
states then one will look for operator valued solutions of the reflection equation so as 
to be able to describe processes where the reflection of a soliton rotates the degenerate 
boundary states among themselves. 

The sinh-Gordon model is just the simplest member of the family of affine Toda 
field theories. For all Toda theories there exist integrable boundary conditions of the 
form p 

r 

a.0(O) = 5^ A,a,e"'-^»/^ (47) 

i=0 

However the boundary parameters Ai were found to be restricted by the requirement of 
integrability to only a small discrete set of values. It would be nice to see if one could 
replace the fixed parameters by the dynamical variables of a mechanical system in 
the way in which we have done for the sinh-Gordon model in this paper. The fixed 
boundary conditions might then arise as the possible stationary points of the boundary 
system. 

Besides the affine Toda field theories, whose integrability is based on the trigonomet- 
ric i?-matrices, there are many integrable field theories related to rational i?-matrices, 
for example the non-linear Schrodinger model and the principal chiral models. Finding 
dynamical solutions of the refiection equation for these rational models is simpler than 
in the trigonometric case and many are already known [O. These could be used to 



construct integrable couplings to boundary mechanical systems for these field theories. 



Acknowledgements: We would like to thank Chris Fewster, Evgueni Sklyanin and 
Anastasia Doikou for very helpful discussions and Vadim Kuznetsov for showing us 
his solutions of the refiection equation. PB is supported by Marie Curie fellowship 
HPMF-CT-1999-00094 and GWD by an EPSRC advanced fellowship. 



References 

[1] Z. Bajnok, L. Palla, G. Takacs, G.Zs. Toth, The spectrum of boundary states in 
sine-Gordon model with integrable boundary conditions; liep-th/010607Cl| . 



14 



R. J. Baxter, "Exactly solvable models in Statistical Mechanics" , Academic Press, 
(1982). 

V. Bazhanov, S. Lukyanov, A. Zamolodchikov, On nonequilihrium states 
in QFT model with boundary interaction, Nucl.Phys. B549 (1999) 529-545; 
tiep-th/98T2091 . 



l.V. Cherednik, Factorizing particles on a half line and root systems, Theor. Math. 
Phys. 61 (1984) 997. 

E. Corrigan, P.E. Dorey, R.H. Rietdijk and R. Sasaki, Ajfine Toda field theory on 
a half line, Phys. Lett. B333 (1994) 83-91; [hep-th/940"4T08| . 



E. Corrigan, Integrable field theory with boundary conditions. Frontiers in Quan- 
tum Field Theory, eds. C.-Z. Zha and K. Wu, (World Scientific 1998), pp 9-32; 
|iep-th/9612138| . 

E. Corrigan and G.W. Delius, Boundary breathers in the sinh-Gordon model, J. 
Phys. A 32 (1999) 8601-8614; |hep-th/9909145i 

E. Corrigan and A. Taormina, Reflection factors and a two-parameter family of 
boundary bound states in the sinh-Gordon model, J. Phys. A 33 (2000) 8739-8754; 
|iep-th/0008237 . 



G.W. Delius, Soliton-preserving boundary conditions in affine Toda field theories, 
Phys. Lett. B444 (1998) 217-223; |hep-th/9809140| . 



L.D. Faddeev and L.A. Takhtajan, Hamiltonian methods in the theory of solitons. 
Springer- Ver lag, Berlin, Germany (1987). 

H. Frahm and N.A. Slavnov, New solutions to the Reflection Equation and the 
projecting method, J. Phys. A 32 (1999) 1547-1555; |cond-mat/9810312| . 

M. Fuentes, A. Lopez, E. Fradkin and E. Moreno, Bosonization rules in 1/2 + 1 
dimensions, Nucl.Phys. B450 (1995) 603-640; |cond-mat/9502076[ 

S. Ghoshal and A.B. Zamolodchikov, Boundary S matrix and boundary state 
in two-dimensional integrable field theory. Int. J. Mod. Phys. 9 (1994) 3841; 
hepth/9306002. 

S. Ghoshal, Bound State Boundary S-matrix of the sine-Gordon Model, Int. J. 
Mod. Phys. A9 (1994) 4801-4810; |hep-th/9310T8^ . 

T. Hollowood, Solitons in affine Toda field theories, Nucl. Phys. B384 (1992) 
523-540. 



15 



V.B. Kuznetsov and A.V. Tsiganov, Separation of variables for the quantum rel- 
ativistic Toda lattices; [hep-th/ 94021 1 1] . 

V.B. Kuznetsov, New Boundary Conditions for Integrable Lattices, J.Phys. A 28 
(1995) 4639-4654; |tiep-th/9503168[ 



V.B. Kuznetsov and A.V. Tsyganov, Quantum relativistic Toda lattices, J. Math. 
Sci. 80 (1996) 1802-1810. 

A. Maclntyre, Integrable boundary conditions for classical sine-Gordon theory, J. 
Phys. A 28 (1995) 1089-1100; |hep-th/9410026i 

P. Mattsson and P. Dorey, Boundary spectrum in the sine-Gordon model with 
Dirichlet boundary conditions; |liep-th/0008071 . 



R.I. Nepomechie, The boundary supersymmetric sine-Gordon model revisited, 
Phys. Lett. B 509 (2001) 183-188; [hep-th/0103029. 



H. Saleur, S. Skorik and N.P. Warner, The boundary sine-Gordon theory classical 
and semi-classical analysis, Nucl. Phys. B441 (1995) 421; |hep-th/9408004 

E.K. Sklyanin, Boundary conditions for integrable equations, Funct. Anal. Appl. 
21 (1987) 164. 

E.K. Sklyanin, Boundary conditions for integrable quantum systems, J. Phys. A 
21 (1988) 2375-2389. 

E.K. Sklyanin, Exact quantization of the sinh-Gordon model, Nucl. Phys. B326 
(1989) 719. 

Yu.B. Suris, Discrete time generalized Toda lattices: complete integrability and 
relation with relativistic Toda lattices, Phys. Lett. A 145 (1990) 113-119. 

V.O. Tarasov, The integrable initial-boundary value problem on a semiline: non- 
linear Schrddinger and sine-Gordon equations. Inverse Problems 7 (1991) 435-449. 

Y.P. Wang, J.H. Dai, Z.N. Hu and F.C Pu, Exact results for a Kondo prob- 
lem in a one- dimensional t-J model, Phys. Rev. Lett. 79 (1997) 1901-1904; 
|cond-mat/970608e . 



A.B. Zamolodchikov and Al.B. Zamolochikov, Factorized S-matrices in two di- 
mensions as the exact solutions of certain relativistic quantum field theory models, 
Ann. Phys. 120 (1979) 253. 

H.-Q. Zhou, X.-Y. Ge, J. Links and M.D. Gould, Graded reflection equation 
algebras and integrable Kondo impurities in the one- dimensional t-J model; 
|cond-mat/9809056| . 



16 



